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We apply the DRA method to the calculation of the four-loop 'QED-type' tadpoles. For arbitrary 
space-time dimensionality T> the results have the form of multiple convergent sums. We use these 
results to obtain the e-expansion of the integrals around T> — 3 and T> = 4. 



I. INTRODUCTION 

The calculation of the high-order radiative corrections have become necessary in many areas of physics, from solid 
state physics to quantum electrodynamics (QED) and quantum chromodynamics (QCD). The radiative corrections 
are expressed in terms of the loop integrals, therefore it is necessary to be able to calculate them. Several powerful 
approaches to this problem have been developed. One of the most successful approach is based on the integration- 
by-parts (IBP) reduction procedure [3, 0. This method allows one to reduce the problem of calculation of arbitrary 
loop integral to that of some finite set of master integrals. The important feature of the IBP reduction is that it 
may also help in the calculation of the master integrals. Namely, using the reduction one can obtain the differential 
and difference @, @] equations for the master integrals. Recently, in Ref. @, the method of calculation based 
on the dimensional recurrence relation Q and analyticity with respect to space-time dimensionality T> (the DRA 
method) was suggested. This method was applied to the calculation of master integrals for several physical problems 
[Tol - [l3j . In the se p apers the DRA method was combined with other methods such as the method of Mellin-Barnes 
representation |l4| and sector decomposition method, implemented in FIESTA (l5| . 

In the present paper we apply the DRA method to the calculation of the 'QED-type' four-loop master integrals 
depicted in FigJlJ These integrals were considered in many papers, see Refs. l9Ul2[|l6l - |25j and references therein. The 
numerical results for their e-expansions around V = 4 were obtained in Ref. [19| using Laporta's difference equation 
method [8( . Some of the integrals are known in analytic form in terms of the hypergeometric function. The integrals 
J6,3 and Jg i have been already investigated using DRA method, see Refs. [t| Q3, HE]. For the maj ority of the 
integrals, several terms of the e-expansion around V = 3 and V = 4 were found in analytical form in Refs. [1911201 l24j . 
However, the complete set of the analytical results for all integrals was not obtained so far. In particular, there is no 
analytical results for the e-expansion of the most complicated non-planar integral Jg ; i around T> — 3. The goal of the 
present paper is twofold. First, we demonstrate some peculiarities of the application of the DRA method appearing 
in the calculation of the integrals considered. Second, we present the complete set of the analytical formulas for the 
e-expansions of the integrals in FigQ] around T> = 3 and V = 4. 

II. THE METHOD OF CALCULATION 

In order to calculate the integrals depicted in FigQ] we use the DRA method [9J , based on the dimensional recurrence 
relation and analytical properties of loop integrals as functions of V. We evaluate integrals in the order determined by 
their complexity level starting from c.l. — 1 and ending with c.l. = 3. The calculation of the integrals J^^fJ)) and 
^7,2 (2?) demands a slight extension of the approach of Ref. |9|], we demonstrate it by presenting here the calculation of 
the integral J-!^{V). Due to the chosen order of calculation, all simpler master integrals of J~j^{V) (the ones obtained 
by contracting some lines of Ji t i{T>)) are already known at this stage. 

Obviously, the integral J~j^{T>) has no infrared divergences for V > 2. Similarly, it has no ultraviolet divergences 
for V < 7/2. Therefore, the integral is a holomorphic function in the stripe {T>\ Re£> £ (2, 7/2)). However, this stripe 
is too narrow to be chosen as a basic stripe of the DRA method. As it was pointed out in Ref. [9] , in this case one can 
pass to a new master integral with dots on the massive lines in order to improve the ultraviolet behavior. A suitable 
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FIG. 1: Four-loop 'QED-type' master integrals considered in this paper. The dashed lines denote massless propagators 1/fc 2 , 
solid lines denote massive propagators l/(k 2 + 1). 




FIG. 2: The integral Jr,i(T>). The dots on the solid lines denote squared propagators. 



choice is the integral J 7 1 (T>) depicted in Fig. [2j The integral 1 (T>) has no infrared divergences for 2? > 2, and no 
ultraviolet divergences for 2? < 9/2, thus, being a holomorphic function in the stripe {2? Re2? e (2,9/2)}. Due to the 
IBP identities, the integral J7,i(2?) can be expressed via Jf 1 (T>) and simpler integrals as 



J 7tl (V) = G(V)+H{V) 



8(225-9) 



3(2?-4) 2 (32?- 10) 



J$AV)+H(p), 



(1) 



where H(T>) is the linear combination of the integrals J4 : i(2?), J^^iV), Js^i'D), J5,3{D), Je t 2(T>), see Appendix. 
Therefore, we can apply the DRA method to the calculation of Jy 1 (2?), and then use Eq. ((T|) to determine Jt^(T>). 
However, we find it convenient to apply the DRA method directly to the integral Jf t x(T>), using Eq. ([lj for the 
determination of the analytical properties of this integral in the basic stripe S — {T>\ KeV s (2, 4]}. The singularities 
of Jt : \(T>) in S are determined by those of two terms in the right hand side of Eq. ([IJ. The singularities of the 
function H(T>) are located at T> = 10/3, 7/2, 4 and totally fixed by the explicit form of the simpler master integrals 
(already calculated at this stage). The only singularities of the first term are the first- and the second-order poles at 
T> = 10/3 and T> = 4, respectively. The principal parts of the Laurent series expansion of the first term around these 
two points are not known. However, by the proper choice of summing factor E(2?) we can achieve that the product 
S(2?)G(2?) is holomorphic in S, see below. 

The dimensional recurrence relation for the integral Jj y \(T>) reads 



J 7) i(£> + 2) 



64(22?-7)(22?-5) 



3(2> - 2) 2 (2? - 1)2?(32? - 8) (32? - 4) 



J 7 .i{V)+R{V) 



(2) 
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where the function R(T>) is given in the Appendix. Choosing the summing factor S(2?) in the form 

2 7 - 2p r(f -p)r(f) 
E(2?) ~ n^r (5 - r(3 - v) ' (3) 

we can rewrite Eq. ^ as 

g(V + 2)=g(V)+r(V), (4) 

where g(D) = £(£>) J 7 ,i(D), r(D) = E(X» + 2)R(V). It follows from Eq. ©, that the function S(D) has first-order 
zeros at 2? = 3, 10/3, second-order zero at T> = 4, and behaves at ImT> — > ±oo as exp{7r|2?|/2}. Therefore, the 
product £(2?)G(£>) is a holomorphic function in S and falls off exponentially at Im£> — > ±oo. Our choice of £(£>) 
provides that the function g(T>) has only known singularities in S, and falls off exponentially at Im£> — > ±oo. In order 
to represent the general solution of Eq. ^ in terms of infinite series, see Ref. [9j, we need to decompose the function 
r(V) as 

r(V) = r+{V)+r_(V + 2), (5) 

where r±(T> ± 2k) decreases faster then 1/fc at fc — >• oo. However, the terms in r(T>), proportional to Je^iD) and 
^5,3(2^), decrease as 1/fc at fc — > 00, and the decomposition (J5]) is not possible. In order to deal with this problem we 
use the following trick. Let us consider the dimensional recurrence relation 

g(V + 2) = g(p) + f(V) , (6) 

for the linear combination 

g(V) = E(X>) (J 7 s(V) + a(V)J 6 , 2 (V) + 0(p)J 6t3 (V)) , (7) 

and try to find the rational functions a(T>), /3(T>) such that the terms proportional to Jq^'D) and J5,3(£>) in r(T>) 
decrease faster than 1/fc at large k. Using the explicit form of r(T>), see Appendix, we choose a(T>) — — 1, and 
conclude that there is no proper choice of /3(X>), and we simply put (3(T>) = 0. After this, the function f(T>) can be 
presented as 

r(V) = r + (V)+f-(V + 2)+f (V), (8) 

where r±(T> ± 2k) decrease faster then 1/fc at fc — >• 00. The function fo(V) corresponds to the large-£> asymptotic of 
the term proportional J^,^{T>) in r(D). Obviously, r (D ± 2k) ~ 1/fc at large fc, so that the sum 

00 

5^fb(2?±2fc) (10) 
diverges. Fortunately, the solution of Eq. (|6]) with r replaced by fo can be written explicitly as 

sin^(7rP/2)^ ( 2" 

where ip( x ) = r'(x)/r(x). Therefore, the general solution of the dimensional recurrence relation (jB]) has the form 



9o(V) = -—27-^-^ ( - ) > (11) 



§(D) ^(z)-^ + (P + 2fc)+^M2?-2fc)+.g (P), (12) 

fc=0 fc=0 

where ui(z) = ui (exp[i7r2?]) is arbitrary periodic function. It follows from Eq. (| 12f) that the analytical properties of 
the function lo(z) are determined by those of g, r±, and go- Namely, u>(z) is a meromorphic function, which has poles 
at z = ±i , ±1, and falls off at \z\ — > 00. Since the principal parts of the Laurent series expansions around these 
points are determined by known integrals only, and not by Jj^x, the function uj{z) can be easily found (for brevity, we 
do not present its explicit form here). Finally, using Eqs. flU, 0, an d (fl2|) we obtain 

J7,x{v) = j 6 , 2 (p) + s- i (p) L(p)-f;f + (p+2fc)+f;f-(p-2fc)- . • (13) 

V ^0 ^0 sm (^p/2) v 2 yy 
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This equation is valid for arbitrary T>, and, in particular, can be used for the numerical calculation of e-expansion of 
Jt t i(T>) around V = 3, 4. In two following sections we present such expansions for all integrals, depicted in FigQ]The 
analytical form of the expansions is obtained from the high-precision numerical results using PSLQ algorithm Ref. 26j , 
as implemented in MPFUN multiple-precision subroutines [27|. The coefficients in the expansions are expressed in 
terms of the following transcendental numbers: 

oo oo oo m / iN TO +fc 

C« = Ef' a » = E^F= Li «( 1 /2), *6= ££^^= 0.98744... (14) 

fe=l fc=l m=lfe=l 

III. EXPANSION AROUND D = 4 
The integrals with c.l.= 



J4,i(4-2e) 

r 4 (-i + e ) = (15) 

J 5 , 3 (4-2e) = 2 1 - 2 T(2-e)r(e^l)r(3e-2) 
r 4 (-l + e) r(e - l)r (2e - |) 

J 5)4 (4 - 2e) _ 3 ■ 2 3 - 4e (e - l)r(2 - e) 2 r (e + |) r(3e - 2)r(4e - 3) 
r 4 '(-l + e) ~ ~~ r(e) 3 r(3e-i) 

J M (4 - 2e) _ r(2 - 3e)r(l - e) 4 r(e) 2 r(3e - l) 2 r(4e - 2) 
r 4 (-l + e) ~~ r(2-2e) 2 r(2-e)r(e-l) 4 r(6e-2) 

The integrals with c.l.— 1 



(16) 



(17) 



(18) 



-V 4 -*) = ? £^ i-4 + ^-^ + ^^ + ^ 2 ln 2 2-^-512a 4 

r 4 (-l + e) (l-3e)(2-3e)(l-2e) | 3 3 V 45 3 3 



/544 4l „ 128 2l 3 „ 1281n 5 2 o „ w „ „,„„ > \ 7 /64tt 6 544 4l 2 „ 

tt 4 ln2 H 7T 2 ln 3 2 h 3072a 5 - 2480C 5 + e 7 1 tt 4 In 2 2 

V 15 3 5 / V 5 5 



+64^ 2 ln 4 2 - _ i 8 432a 6 - 7680s 6 + + (e 8 ) } , ( 10 ) 



■/ 5 ,2(4 - 2e) (1-e) 3 I_fi + ™ 344e 2 3584e 5 C 3 6 /87047T 4 

T 4 (-l + e) (l-4e)(3-4e)(l-3e)(2-3e)(l-2e) \ 3 3 6 \ 45 

2048 2n 2048 ln 4 2 \ 7 f 34816 4l „ 8192 2l , 8192 ln 5 2 

H 7r 2 ln 2 2 16384a 4 + e 7 tt 4 ln2 H 7r 2 ln 3 2 

3 3 / \ 15 3 5 

\ o/13312tt 6 69632 , 2 , , 16384 ln 6 2 

+196608a 5 - 174592C 5 - e 8 + tt 4 ln 2 2 + 8192tt 2 ln 4 2 

/ \ 9 5 5 

-2359296a 6 - 1081344s 6 + 1266 ^ 88 Cl \ + Q ^ I ; (20) 



J 62 (4-2e) 2 4e 2e 2 ,/44 I6C3 \ 4/ 200C 3 4tt 4 \ ./1928 1192C 3 326tt 4 

-el— — - e 116 h — — — e 1 



r 4 (-l + e) 3 3 3 V 3 3 / V 315/ V 3 3 45 

64 2l 2 „ 641n 4 2 „„ \ fi /9328 5864C 3 2126tt 4 448 2l 2 „ 

_ 3584oi - ™* ,„ 2 - ,„3 2 + _ 4096 „ 5 + 2784C5 + £ 

o 45 9 15 21 

+<*'). 



r 4 (-l + e) 4 2 V 2 / V 4 8/ V 2 5 2 

The integrals with c.l.= 2 



J M (4-2e) _ 3 , 7e , 9e 2 ^39 , \ _ ^ / ^ _ + ^ _ ^ / ^ _ , 54T. 4 



T 4 (-l + e) 2 2 2 V 2 7 V 2 0/ V 60 

- e 6 ( 6336 - 4851C 3 + — ^— + 2 40tt 2 ln 2 2 

272 1921n 5 2 177r 6 
-240 ln 4 2 - 5760a 4 - — tt 4 ln2 - 64tt 2 ln 3 2 + 4608a 5 + 3531C 5 + -^j- + 498^ 

2i 1 

+0(e 7 ), 



r 4 (-l + e) 6 6 V 3 J V 3 3 60 y V 3 3 6 

5 /602 38C 3 85tt 4 16 2 , 2o 161n 4 2 nn irj> 44tt 6 „ 

-e 5 ^ + - r— - + — tt 2 ln 2 2 I28a 4 - 154C 5 + —^r + 128C? 

V 3 3 36 3 3 s 189 S3 



J 7 , 2 (4-2e) _ l 5e ^ /ll Cs \_ £ s f 44 l3 Cs , tt 4 ^ 4 fl66 29C 3 , 5tt 4 43( 5 



r 4 (-l + e) 66 V 32 / V 3 6 120/ V 3 6 24 

_ 5 /602 197C 3 41tt 4 23lC 5 17tt 6 105C 3 2 \ ^, 6 , 
6 \ 3 6 120 2 189 2 / 16 j ' 



The integral with c.l.— 3 



124 4l ,„ 80 2l 4o 321n 6 2 2859C?\ ^. fiN 
— — — 7r ln 2 2 - — tt 2 ln 4 2 + + 7680a 6 + 4032s 6 — ^ + (9(e 6 ) 
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The above expansions were considered in Refs. (13, Il6l - [24j . Our result for Jg j2 is in full agreement with that of Ref . (2l| . 
The analytical form of the expansion of all integrals, except the most complicated integral Jg \, up to the terms with 
maximal transcendentality weight equal to 5 was presented in Ref. (l9j |. The numerical form of the expansion of Jg^ 
was calculated in the same paper, however, the precision of this calculation was not sufficient for the application of the 
PSLQ algorithm. The e° term of Jg i in analytical form was calculated in Refs. (HI, In Ref. [24[ the e-expansion 
around T> = 4 of some integrals was presented up to the terms with maximal transcendentality weight equal to 8. 
Our exact expressions can be immediately used for the extraction of even higher terms of e-expansion, but we assume 
that the practical significance of these terms is questionable. 



IV. EXPANSION AROUND D = 3 



The integrals with c.l.= 1 



J 5 i(3-2e) (1 - 2e) 2 J 1 Jln /2tt 2 Jl2n \ 2 /4 2l „ 81n 3 2 on \ ,/44tt 4 
- 4 ln2 + e I — - — h 4 In 2 I -e —tt ln2 H h 38C 3 + 6 



r 4 (-l/2 + e) l-6e 1 e \ 3 J \3 3 "V V 45 

16 9 2 4 A 4/^88 , 32 , o 161n 5 2 

7T 2 In 2 + 8 In 2 + 160a 4 + 216 ln2 (3 - e — tt 4 ln2 tt 2 ln 3 2 + 320a 5 

3 / \45 9 5 

+36^ 2 C 3 + 216 ln 2 2 C 3 + 1M5{ 5 ) + e 5 ( + ^tt 4 ln 2 2 - ±?7r 2 ln 4 2 + + 640a 6 



V 27 45 9 15 

1568s 6 + 727r 2 ln2C 3 + 1441n 3 2C 3 + 1614Cf + 59281n2C 5 ) + 0(e 6 ) I , (28) 



« 3 -^ - (1 - 2£)3 .{--81u2 +e ^ + 161n 2 2V e 2 f-- 2 ln2 + ^ + 108 C3 



r 4 (-l/2 + e) (1 -6e)(l -4e) |4e 



3 



3 /316vr 4 „„ 2l 2 „ 801n 4 2 „„ n „ \ 4 /1264 4l „ 64 2l , r 

-e 3 h 16tt 2 ln 2 2 + h 128a 4 + 544 ln2 £ 3 - e tt ln2 H tt 2 ln 3 ^ 

\ 45 3 / \ 45 3 

641n 5 2 544tt 2 C 3 2 \ ,/21928tt 6 2528 , 2 

512a 5 + + 1088 ln 2 2 ( 3 + 3212C 5 + e n „ K + In 2 

3 3 /V 945 45 



64 2l 4 „ 1281n b 2 „„„ 2176 2l „ „ 4352 

-— tt 2 ln 4 2 + + 2048a 6 - 256s 6 + -^-tt 2 ln2 C 3 + -r- 

o 9 00 



+13344 ln2C 5 J + 0(e b )j , (29) 
J 6 , 2 (3-2e) _ n ,3/^ / 1 2 7C 3 \ , /897T 4 1 2l 3 ln 4 2 \ ,/ 89 4 . 



(1 - 2e) d { — 7r 2 ln2 + — +e tt 2 In 2 + — + 8a 4 -e tt 4 ln2 



r 4 (-l/2 + e) v ' 1 32e \ 8 4 / V 144 12 3 ) \ 360 

,1 2, 3 9 41n 5 2 127^ 2 C 3 403^ 3 ^ 13159^ 6 1 4 2 1 2 4 , 81n 6 2 
+ -7r m 2 h 32a5 H 1 + e 1 7r In 2 H — 7r In 2 1 



9 15 48 16 y V 60480 20 3 45 

327r 2 a 4 5 „ 253C?\ , ,,1 

+— 3-^ + 128a 6 - 52s 6 - -tt 2 ln2 C 3 + -^j + 0(e 4 ) | , (30) 

J 63 (3-2e) . tt 2 / 3 2l „ 21C 3 \ / 23tt 4 3 2l 2 „ 31n 4 2 n „ \ 2 /69 4l „ 

r 4 (-i/2 + e) = (1 - 2e %iy in2+ T) + {-W + f k2 + — + H- e tio /ln2 

-3. 2 ln 3 2 - + 432 fl5 + ^ - ^ + e 3 ( - ^ - ^ 4 ln 2 2 + W ln 4 2 

5 16 16 y V 448 20 

+ 36^ + ^ +1 + 56M^\ +0(£ 4 } \ (31) 

5 4 8 / 



7 



in 2 ( ir 2 11( 3 \ ,/ 27tt 2 2 , „ 51C 3 271tt 4 \ ,/907tt 2 lr 2l „ 

- ^ + ^ ln 2 2 + ^ + 25 C5 ) + , 4 + 129. 2 In2 + ^ 

16 2880 6 s / V 96 16 

2159tt 4 n „ 2l 2 „ 14 4 , „ 8 2l o 37tt 2 C 3 100> 10279tt 6 4l ,„ 

^ 307r ln 2 + ^r 71 " ln2 + 77 71 " ln 2 tt^ ~ 183 Cs + + In 2 

320 3 3 3 12960 

-7T 2 ln 4 2 - 24^ 2 a 4 - 21tt 2 In2 C 3 - + 0(e 5 ) j . (32) 

The integrals with c.l.— 2 



J 6 i(3-2e) , N ,J 7T 2 / 1 2, 21C 3 \ / 5tt 4 2 , , 51n 4 2 63, \ 

1%-1/2 + e) = (1 - 2£) \327- ln2+ ^J +e (-l^- 7r 2 ln 22 + -^ + 30 fl4 + T ln2 C3 j 

£ 2 f -7T 4 ln2 - V ln 3 2 +1^ + 108 ln2 a 4 + 228« 5 + ^ + ln 2 2 C 3 - ^ 
V12 3 5 5 16 4 Sd 8 

+e 3 f- 5 !^ - T^^ 4 ln22 + ln42 + ^ ^ + 18ir 2 a 4 + 108 ln 2 2 a 4 + 648 ln2 a 5 + 1560a 6 
\ 6720 12 6 15 

+336 S6 + 9 -n 2 ln2 C 3 + y ln 3 2 ( 3 - ^ - ^ ln2 C 5 ) + 0(e 4 ) j , (33) 

J 7 i(3-2e) (l-2e) 3 f/7r 2 ln 2 2 \ /3 2n „ , 3 „ i( ^ 23 2, 2„ ln 4 2 „„ 

' 1 e -7T 2 ln2 + ln 3 2 - 4C 3 + e 2 — - — it 2 ln 2 2 + — - + 30a 4 



(l-2e) 3 . 


if- 2 


4e+ 1 




-|ln2 C3 


h 


2103C 5 \ 


+ £ 4 I 


16 ) 



r 4 (-l/2 + e) 4e + l 1 V 24 2 J \4 ™J \ 144 12 12 



361 4 , n 7 2l 3 „ 121n 5 2 „ o , „ „ o 13tt 2 C 3 213, 2n , 

— 7T 4 ln2 + -7T 2 ln 3 2 + + 28 ln2 a 4 - 28a 5 + — ln 2 2 £3 

180 9 5 6 4 

/9361tt 6 79 4 , 2 „ 1 4 „ 1331n 6 2 „ 2 ™, 2„ 

7T 4 ln 2 2 + — 7T 2 ln 4 2 27r 2 a 4 - 92 ln 2 2 a 4 - 184 ln2 a 5 

V 11340 36 18 45 

+24a 6 - 278s 6 -^ 2 ln2 C 3 - ^ ln 3 2 C 3 + - ™ ln2 (5) + 0(e 5 ) \ , (34) 



J 72 (3-2e) . „ N3 [3C 3 / 7T 4 1 2l o ln 4 2 o „, „ _ \ 2 /8 2 , 3 „ 161n 5 2 
1^-1/2 + e) = (1 ~ 2£) (rVF3' 2 ln 22 + - r+ 8 fl4 + 71n2 C3 j +£ 2 l^ 2 - — 

-32 ln2 a 4 - 32a 5 - ^ - 14 ln 2 2 ( 3 + ^ + , 3 f - - V m 2 2 - H. 2 ln 4 2 

5 12 S3 16 J V 45360 9 9 

161n 6 2 87r 2 a 4 2 „ _ o , „ „„ o nA 7 2l 56 , ,„ , 181(1 

H 1 -— - + 64 ln 2 2 a 4 + 128 ln2 a 5 + 128a 6 - 84s 6 + -tt 2 ln2 ( 3 + — ln 3 2 C 3 + 



3 3 ' 3 ' 2 

- ln9 i^r I 



+ ^iln2C 5 ) +0(6 4 )^. (35, 
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The integral with c.l.— 3 



J 9il (3-2e) ( 3 31n2 3tt 2 ( 91n 2 2 C3 ^\ ( 417 309 ln2 i 37tt 2 331n 2 2 9 ., 



r 4 (-l/2 + e) \32 128 256 64 64/ V 256 256 512 128 64 



/ 417 


309 ln2 


V 256 


+ 256 


ln 4 2 


21 + 1] 


192 + 


8 64 



tT ln2 

91n 3 2 9C 3 43tt 4 1 2l 2 „ ln 4 2 a 4 7, „ \ 2 /8313 66871n2 

-tt ln 2 1 " 1 



32 8 23040 192 192 8 64 "J \ 512 512 

223tt 2 1653 ln 2 2 29 2 351n 3 2 77( 3 517tt 4 21 2 2 631n 4 2 
-7r m2 H — — — — — + — — -7T In 2 1 



1024 256 128 ' 64 256 5120 128 128 

63a 4 1071 1 2 3 9 ln 5 2 1 a 5 11tt 2 C 3 7 2 251Cs 

ln2 C3 H 7T In 2 m2 a 4 In 2 c 3 H 

16 128 s 144 120 4 4 128 64 S3 256 

,2, 



,/ 131151 80349 ln2 2461tt 2 27903 ln 2 2 2295 2l 1447 ln 3 2 18557C 3 

-e 3 7r 2 ln2 1 — 

V 1024 1024 2048 512 256 128 512 

834197T 4 1435 2l 2 „ 3199 ln 4 2 2701a 4 13249, „ „ 17 4l „ 29 2l 3 „ 

7T 2 ln 2 2 H 1 H ln2 C 3 H tt 4 ln2 tt 2 ln 3 2 

92160 768 768 32 256 s 40 32 

391n^2 405 423as 699^ 4311 2 63189C 5 1445^ 6 17 4 2 

+ ^^ + — ln2a4 + ^ + ^5^ + i2y ln2C3 -^T2-- ^Tey-^ ln2 

j_ 3 2 4 ln 6 2 17A 1 2 1 a 6 23s 6 119 2 7 3 

H 7r In 2 ln 2 a 4 ln2 05 n ln2 c 3 In 2 C 3 

32 144 8 4 2 2 8 64 s 96 

Some of the above expansions were considered in Refs. Hp . I20L [24l |2^ |. The e-expansions for the integrals and 
Jq.3 are in agreement with those obtained in Refs. [13, l28| up to the terms considered in these papers. In Ref. [24| 
some higher terms of the e-expansions of the integrals Js^, Js,2, J&,2, ^6, 3l ^7,1 were presented. However, we observed 
several inconsistencies of the results of Ref. [24[ with our results. In particular, the e 4 terms in Eqs.(5.10), (5.22), and 
(5.28) of Ref. Q seem to be incorrect. 



V. CONCLUSION 



In the present paper we apply the DRA method to the calculation of the four-loop 'QED-type' tadpole master 
integrals. The results obtained are valid for arbitrary T>, and have form of the convergent multiple sums. For brevity, 
these results are not presented here, and are available from the authors upon request. We have presented the e- 
expansions of the integrals around T> — 3, 4. The highest transcendentality weight of the expansions (equal to 6), 
was chosen rather arbitrarily and should be sufficient for physical applications. Higher terms of e-expansion for all 
considered integrals can be easily obtained from our exact in T> expressions for the integrals. 
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VI. APPENDIX 



= (P-3) (52? 3 - 592? 2 + 2302? - 297) (X> - 3)(X> - 2)(32? - 8) (2? 2 - 122? + 30) 

1 j 3(2? - 4) 3 (32? - 10) 6 ' 21 ' 24(2? - 4) 3 (22? -7)(32?- 10) 5,31 ' 

(22? - 5) (325 - 11) (32? - 8) (42? 2 - 292? + 54) 

96(2? -4)3(2? -3) (32?- 10) J W V ) 
(2? - 2)(32? - 8) (182? 3 - 2152? 2 + 8452? - 1098) 
+ 96(2? -4)3(2? -3) (32?- 10) ■ / m( i> ) 

(2?-2)3(42?3 + 152? 2 - 2292? + 450) 

192(2? - 4)3(2? -3) 2 (32?- 10) j4 > 1 ^' ^ 

J2(X>) = A( 6 ' 2 )(2?)J 6i2 (X?) + A( 5 ' 3 )(2?)J 5 , 3 (I?) + ^ 5 ' 2 )(2?)J 5 , 2 (2?) +A( 5 ' 1 )(2?)J 5;1 (2?) +A( 4 ' 1 )(2?)J 4; i(2?) , (38) 
where the coefficients A^ % '^ (2?) have the form: 

(6)2) 16(2? - 3) (372? 4 - 2692? 3 + 6892? 2 - 7182? + 246) 

^ = 3(2? - 2) 2 (2? - 1) 2 2?(22? - 5X22? - 3) (32? - 8) (32? - 4) ' 

(5 3) _ 8 (5672? 7 - 83702? 6 + 524452? 5 - 1806392? 4 + 3690212? 3 - 4466962? 2 + 2964002? - 83088) 

^ ' ~ 9(2? - 3)(2? - 2) 2 (2? - 1) 2 2?(22? - 5) (22? - 3) (32? - 8) (32? - 4) 2 ' 

(5)2) 8 (2162? 4 - 17922? 3 + 55152? 2 - 74792? + 3780) 

' ~9(2? - 3) (2? - 2) (2? - 1) 2 2?(22? - 3) (32? - 8) (32? - 4) 2 ' 
(5 1) 16 (722? 4 - 5762? 3 + 16912? 2 - 21712? + 1044) 

(2?) ~ 9(2? - 3)(2? - 2)(2? - 1) 2 2?(32? - 8)(32? - 4) 2 ' 
(41) 4 (92? 6 - 4232? 5 + 35272? 4 - 125602? 3 + 224492? 2 - 198542? + 6912) 

A ' (P) = 9(P - 3) 2 (2? - 2)(2? - 1) 2 2?(22? - 3)(32? - 8)(32? - 4) 2 ' ^ 39 ^ 

f(2?) = S(P + 2)|(A (6 ' 2) (2?) -a(P)yl (7 ' 1) (P) + a(2? + 2)B (6 ' 2) (2?)) J 6>2 (2?) 

+ (a^{V) + a(2? + 2)B (5 ' 3 \V) - /3(2?)A^ 1 )(2?) + /3(2? + 2)C( 5 < 3 )(2?)) J 5:3 (2?) 
+A( 5 ' 2 )(2?)J5, 2 (2?)+A( 5 ' 1 )(2?)J5, 1 (2?)+A^ 1 )(2?)J4,i(2?)l , (40) 



where 



A{ 7.i) (v) = 64(22?-7)(22?-5) 

1 ; 3(2?-2) 2 (2?-l)2?(32?-8)(32?-4) ' 

B (e,2 )(v) = 16(2?-3)(2?-2) 

k ' (2?- l) 3 2?(22?-5)(22?-3) ' 

(5j3) 8 (372? 4 - 2862? 3 + 8112? 2 - 9962? + 446) 

^ 3(2? - 2)(2? - 1) 3 2?(22? - 5) (22? - 3) (32? - 4) ' 

(2?) = 64(22?-5)(22?-3) 

V ' 3(2?-2)(2?- l)2? 2 (32?-4)(32?-2) V ' 
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